Lineshape distortion in a nonlinear auto-oscillator near generation threshold: 
Application to spin-torque nano-oscillators 
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The lineshape in an auto-oscillator with a large nonlinear frequency shift in the presence of thermal 
noise is calculated. Near the generation threshold, this lineshape becomes strongly non-Lorentzian, 
broadened, and asymmetric. A Lorentzian lineshape is recovered far below and far above threshold, 
which suggests that lineshape distortions provide a signature of the generation threshold. The 
theory developed adequately describes the observed behavior of a strongly nonlinear spin-torque 
nano-oscillator. 
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Thermal noise plays a key role in the dynamics of auto- 
oscillatory systems, especially for the nano-sized systems 
in which the characteristic energy of auto-oscillation is 
comparable with the thermal noise energy ksT. One 
important consequence of noise is the broadening of the 
power spectrum. It is well-established, for example, that 
white Gaussian noise leads to a broadened Lorentzian 
shape of the power spectrum both far below and far above 
the generation threshold in conventional auto-oscillators 
(for which the frequency is weakly dependent on the am- 
plitude) 1]. In addition, the presence of noise also blurs 
the threshold between the damped and auto-oscillatory 
states. 

In the immediate vicinity of the auto-oscillation thresh- 
old, a non-Lorentzian spectrum appears: the lineshape 
can be described by a sum of several Lorentzian profiles 
corresponding to several simultaneously excited statis- 
tical modes. These modes have the same central fre- 
quency for a conventional auto-oscillator and, therefore, 
the contributions of higher-order modes do not lead to 
large qualitative changes in the overall spectral line pro- 
file fH. In contrast, the central frequencies of the par- 
tial Lorentzian profiles are shifted in frequency by dif- 
ferent amounts in a nonlinear auto-oscillator in which 
the frequency is strongly dependent on the amplitude. 
In the vicinity of the auto-oscillation threshold at which 
these shifts are of the order of the oscillation linewidth, 
the frequency nonlinearity leads to important qualitative 
changes in the lineshape: the spectral line is asymmetric 
and significantly broadened. 

An example of an auto-oscillating system with strong 
frequency nonlinearity is the spin-torque nano-oscillator 
(STNO), which is of considerable interest at present. 
This oscillator consists of two ("free" and "fixed") fer- 
romagnetic layers separated by a nonmagnetic (metallic 
or dielectric) spacer. When driven by a spin-polarized 
direct current, an effective negative damping results in 
the "free" layer due to the spin-transfer effect [1, 



For a sufficiently large direct current, the induced nega- 
tive damping overcomes the natural positive damping of 
the "free" layer, and auto-oscillations of magnetization of 
this layer result 0, 0, @] ■ These and related phenomena 
are subjects of intensive research at present due to their 
possible applications in microwave nano-electronics. 

In this Letter, we present a theory of lineshape distor- 
tion in auto-oscillators with a strong nonlinear frequency 
shift in the presence of thermal noise. We show that line 
broadening and line asymmetry are key signatures of the 
threshold region, and we argue that such features provide 
a simple way to determine the threshold in experiment. 
This point is especially significant for STNOs in which 
the influence of noise on the device performance is of a 
qualitative importance. The theory is compared to re- 
cent experiments on STNO nanopillars for which good 
qualitative agreement is found. 

Our theory is based on a classical Hamiltonian for- 
malism for spin-waves in magnetic multilayers with the 
spin-transfer effect 0, It is assumed that only one 
spin-wave mode - the mode having the lowest relaxation 
rate - is excited at the auto-oscillation threshold, and 
dominates the oscillation dynamics both below and just 
above this threshold. Using this assumption and follow- 
ing Refs.0,i we can reduce the Landau-Lifshitz equation 
with Gilbert damping, which describes the magnetization 
dynamics of the STNO "free" layer, to a stochastic non- 
linear oscillator equation for the complex arnplitude c of 
the single excited spin-wave mode (|c| < 1) [lO|, 



dc 
dt 



Here, 
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+ iuj{p)c + r+{p)c - r_(/,p)c = fit). (1) 



in the dimcnsionless oscillation power 



which determines the polar angle 6 of the magnetiza- 
tion precession in the "free" layer 6 = arccos(Af2/Mo) = 
arccos(l — 2|cp), where Mq is the length of the magneti- 
zation vector and is the projection of this vector on 
the direction of stationary equilibrium magnetization z 



2 



(see Ref. for details). The oscillation frequency uj{p) is 
given by uj{p) ~ + Np, where luq is the (linear) mode 
frequency, and the coefficient N describes a nonlinear fre- 
quency shift. Typically for STNOs, |A^| > ujq, which is 
the key feature distinguishing STNO from other conven- 
tional "quasi- linear" auto-oscillators for which iV ~ 0. 
The second term in Eq. ([T]), r+(p) ~ ro(l -I- Qp), repre- 
sents natural positive spin- wave damping, where Tq is the 
linear relaxation rate and Q is the phenomenological di- 
mensionless parameter characterizing nonlinearity of the 
damping yjj. The third term, r_(/,p) = al{l — p), 
describes current-induced nonlinear damping (which is 
negative for small oscillation amplitudes), where / is the 
bias current, and a is the spin-polarization efficiency de- 
fined by Eq. (5) in Ref. 0. 

The function f{t) in (IT|) is a stochastic term that de- 
scribes the influence of the thermal fluctuations. We 
take f{t) to represent a white Gaussian process with zero 
mean and spectral properties of the form {f{t)f*{t')) = 
2TQrj5{t — t'). The strength of the stochastic term is 
given by the product of the linear relaxation rate To 
in the "free" layer and the thermal noise energy r/ = 
ksT / {\u)q)^ where A is the constant that relates dimen- 
sionless oscillation power p = \c\'^ to the energy of the 
system E{p) = XujQp. In our case A = VMq/^, where 
V is the volume of magnetic material involved in auto- 
oscillations, Mo is the saturation magnetization, and 7 
is the gyromagnetic ratio. This choice of noise ampli- 
tude T] ensures that the system, in the absence of spin- 
transfer, relaxes towards the equilibrium thermal noise 
level {E{c)) = ksT at the rate Tq ^,Kd\. 

It should be stressed that the nonlinear oscillator equa- 
tion ([T]) is rather general, and applies not only to spin- 
transfer-induced dynamics in magnetic multilayers, but 
also to the dynamics of any nonlinear auto-oscillating 
system operating in a single-mode regime. 

In experiments involving STNOs the magnetization os- 
cillations driven by spin-transfer torque are characterized 
by the spectral density of the associated voltage oscilla- 
tions. It has been shown elsewhere that the autocorrela- 
tion function of experimentally measured voltage oscilla- 
tions v{t) is directly proportional to the autocorrelation 
function of the corresponding spin-wave oscillations to 
lowest order [9], i.e. {v{t)v{0)) oc Re[(c(i)c*(0))]. In what 
follows, we demonstrate how to calculate the correlation 
fimction K{t) = (c(t)c*(0)) or the power spectrum. 
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of spin- wave oscillations resulting from the stochastic dy- 
namics described by Eq. ([!]) at any magnitude of the 
driving bias current /. 

It is possible to obtain a linearized form of Eq. ^ far 
above and far below threshold, from which the spin- wave 
correlation functions can be evaluated exactly. In the 
subcritical regime, C = ///th <^ 1 (where /th = Fo/cr 



is the threshold current at which the auto-oscillations 
start), all nonlinearities can be ignored and Eq. ([T]) leads 
to a Lorentzian power spectrum S{uj) with linewidth 



Acj = Fn — al. 



(3) 



This linear dependence has recently been observed in 
spin- valve nanopillars 12 1 and tunnel junctions [l^ in 
the low-current regime. Far above the threshold [C, > 1), 
the lineshape is also Lorentzian with a linewidth governed 
by phase noise [1, [3| , 



Au; = Fo (kBT/Eo) [l 



(4) 



where Eo = Xujqp = Ati>o(C ^ 1)/(C + Q) is the av- 
erage energy of the stable auto-oscillation in STNO, 
N = duj{p)/dp is the nonlinear frequency shift coeffi- 
cient, and Feff — T^{p)/dp — dV ^{p) / dp is the effective 
nonlinear damping (see Ref. 14 for details). 

Near threshold, the coupling between amplitude and 
phase fluctuations leads to a non-Lorentzian power spec- 
trum. To demonstrate this, we use the polar coordi- 
nates c{t) — r(t) exp[i(f){t) — iwoi] and linearize ([T]) about 
the steady-state trajectory, ro — \/_{C — 1)/(C + Q)j just 
above threshold, r(t) = ro + a(t) [lj|, 

d(<) + 2reffr2a(i) = Re[/(t)], (5) 
m = -Im[/>)] - 2Arroa(0, (6) 

where /(i) — J(i)e"^f t-i0(t) represents a white Gaussian 
noise with the same statistical properties as j(t). While 
the amplitude fluctuations a(t) are subjected to a white- 
noise forcing only, a(t) itself becomes a colored-noise 
source for the phase fluctuations (/)(t), since {a{t)a{t')) ~ 
exp ( — 2Foffro | i — t' | ) . The stochastic process of the phase 
variable (f)(t) is therefore driven by a colored-noise source 
a{t) and a white noise source f{t), which clearly leads to 
a non-Lorentzian power spectrum if we take the STNO 
to be a simple phase oscillator. 

This linearization scheme is not accurate in the thresh- 
old region because the steady state amplitude ro is com- 
parable to the magnitude of the amplitude fluctuations. 
However, instead of solving the full nonlinear stochastic 
problem ([1]), it is more tractable to solve the correspond- 
ing linear Fokker-Planck equation, which describes the 
time evolution of the probability density function (PDF) 
V{t, c; Co). A similar approach has been used to study the 
problems of thermally-assisted magnetization reversal in 
the presence of spin-transfer torque [3] and telegraph 
noise in bi-stable current-induced precessional states [l6| . 
The PDF V{t, c; co) gives the probability of the auto- 
oscillator to have the complex amplitude c at the time t, 
given the amplitude at the zero moment of time {t = 0) 
was Co. From this non-stationary PDF V{t, c; cq), the 
correlation function K{t) for < > can be obtained as 
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d'^co'Pit,c;co)Vo{co)ccQ, (7) 
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where Po(c) is the stationary PDF, i.e. 'Po(c) = 
Vmit^ac P{t, c; cq). For t < one can use relation 
K(t)^K*{-t). 

The Unear Fokker-Planck equation can be obtained 
from ([T]) in a standard way 17[, which in polar coor- 
dinates is given by 



^ = 0.5 



^= 1.04 



^=1.5 
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r dr 



{[r+(r2) -r_(/,r2)] r^V}+Nr 
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j,2 g^2 



= LV. (8) 



The stationary solution Vq of Eq. ([8]) does not depend on 
the time t and angle </>, i.e. Vq — Vo{r), and can be found 
explicitly in a general case [lO|. Eq. ([8|) is a linear equa- 
tion that does not depend on time t and phase </> explic- 
itly. Therefore, its solution c;co) = V(t,r,4>;rQ,(j)Q) 
can be written as a series 



V{t,r, (jj; ro,(/)o) 



(9) 

where Pn.ii{r) are normalized eigensolutions of the prob- 
lem L„P„,^ = -^n,^iPn,^i, with L„ = L{dr,-in). 
The eigenvalues A„^^ and coefficients a„_^(ro,0o) sat- 
isfy the initial conditions for the PDF: 7^(0, r, 0; ro, 0o) — 
{\/r)5{r — ro)6{(j) — (j)Q). In the expansion ([9]), the index 
n describes the 0-dependence ^ exp(— m0) of the partic- 
ular statistical mode P„,^, while the index fi enumerates 
different modes having the same 0-dependence. 

With the normalized solutions Qn.p(r) of the adjoint 
eigenproblcm, i+Q,i,^ = ~^n,i^iQn,p., and using Eqs. ([7]) 
and ([9]) with Eq. ([2]), a series expansion for the power 
spectrum S{uj) can be obtained in the form. 



Fi,^Re(Ai,^) 



Im(Ai,^)]2 -I- [Re(Ai,^)]2 



(10) 



whereby each contribution is weighted by the factor 
Fi.^ = !^ dr dr' {rr')P,^^{r)Vn{r)Ql^{r'). Each 
term in (|10p is associated with one statistical eigenmode 
(n, /x) . Only the n = 1 modes give contributions to 
the second-order correlator K{t) and power spectrum 
S{uj). Statistical modes with other phase indices n would 
be necessary to calculate the higher-order correlators, 
such as {[c* {t)c{Q)Y) . The partial lineshape associated 
with each statistical mode is Lorentzian, with the cen- 
tral frequency and linewidth given by the imaginary and 
real parts of the eigenvalues Ai_^, respectively. The re- 
sult in pUj) . therefore, represents an expansion of the 
non-Lorentzian spectrum in terms of partial Lorentzians, 
whereby each mode /i > 1 represents an additional re- 
laxation process (and time scale) due to the additional 
colored-noise source. 

A comparison between the lineshape distortion of a 
nonlinear {N = N/Tq = -20) and linear (TV = 0) oscilla- 
tors is presented in Fig. [T] The noise level considered is 
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FIG. 1: (Color online) Power spectrum, S{uj), for r; = 10"'' 
and Q = 2. The dominant fj, modes contributing to the overall 
spectral line(thick black curve) are shown for two values of fre- 
quency nonlinearity iV at different supercriticality C, = I / hi-i- 



rj = 10~'^, which corresponds approximately to T = 300 
K for the nanopillar system in [l8| and T = 400 K for 
the nanocontact system in |^ with an in-plane applied 
field of 0.6 T. Far below and above threshold, the spec- 
tral lines are largely Lorentzian and are well described 
by a single ^ — 1 mode for which the weighting factor 
Fi^i is at least two orders of magnitude larger than the 
others, as one can clearly see in Fig. [1] A significant de- 
formation of the spectral line occurs near and just above 
threshold, as discussed earlier, which appears as large 
contributions from the higher-order eigenmodes in (jlOp . 
For the nonlinear oscillator, there is an additional line- 
shape asymmetry due to the different central frequencies 
of the partial Lorentzians. As each Lorentzian represents 
a stochastic process related to a different amplitude pro- 
cess (resulting in a different relaxation rate), there is an 
associated frequency shift due to the nonlinearity that is 
different for each process. A detailed description of these 
processes will be presented elsewhere. 

Salient features of such lineshape asymmetry have been 

In Fig. m the ex- 
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observed in a recent experiment 
perimental spectral lines are fitted and compared with 
the theoretical lineshapes predicted by the Fokker-Planck 
theory. Based on the linewidth and power variations 
of the experimental data we estimated the experimen- 
tal critical current to be /th = 5.2 mA, which allows us 
to compare the theoretical and experimental lineshapes 
at the same supercriticality C,. Below threshold (/ ~ 5.0 
mA, C « 0.96) the experimental data are well described 
by a single Lorentzian curve (Fig. [5^) , which is consistent 
with the theoretical curve (Fig. [2^) . Just above thresh- 
old (I = 5.4 mA, C, ~ 1-04), a slight asymmetry appears 
in the spectral line whereby a second Lorentzian peak 
is required for a reasonable fit. The presence of a sec- 
ond mode is also seen in the theoretical curve, with a 
good qualitative agreement in the relative positions and 
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FIG. 2: (Color online) Comparison of theoretical and exper- 
imental spectral lines near threshold (Jth ~ 5.2 mA). Thick 
solid lines in (a)-(c) represent spectral lines as predicted by 
Eq. (|10|l . with the first dominant fi modes shown as shaded 
curves. Circles in (d)-(f) are experimental data at T = 225 K 
taken from Ref. 118, with thick solid lines representing a dou- 
ble Lorentzian fit. Shaded curves are individual Lorentzian 
profiles used in the fits. 
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cal regimes. Good agreement between the numerically 
and analytically calculated linewidths is found far below 
and far above the threshold in all cases, as demonstrated 
by the dashed lines representing Eqs. Q and re- 
spectively. Experimentally measured variations of the 
STNO linewidth with supercriticality are in good quali- 
tative agreement with the developed theory (see Ref. To'). 
The variation of the oscillation frequency (Fig.jSb) is rel- 
atively slow below the threshold and exhibits a fast quasi- 
linear decrease with current above the threshold due to 
the nonlinear frequency shift (determined by the coeffi- 
cient N). Only the frequency redshifts {N < 0) have 
been presented here for the sake of brevity; analogous 
behavior is predicted for the frequency blueshifts. 

In summary, we have presented a theory of lineshape 
distortion for auto-oscillators with strong frequency non- 
linearities. The Lorentzian lineshape of the power spec- 
trum becomes asymmetric in the threshold region. Good 
qualitative agreement with recent experimental data on 
the spin-transfer nano-oscillators is demonstrated. 
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FIG. 3: (Color online) (a) Scaled linewidth and (b) frequency 
shift, SlOc — uJc — 1^0, as a function of supercriticality from 
Lorentizian fits to spectra computed for different iV = N/To 
for Q = 2, r]o = 10"'^. Dashed lines in (a) correspond to 
Eqs. © and Q. 



amplitudes of the two modes (Fig. [SJj). At a moderate 
above-threshold current (7 = 6.0 mA, ( « 1.15), the line- 
shape distortion becomes more pronounced (Fig.[2f), and 
this asymmetry is qualitatively reproduced by the theory 
(Fig.H:). 

Another way of quantifying the lineshape distortion is 
provided by single Lorentizian fits to the spectra which 
yield the central frequency Wc and the linewidth Alo. The 
results of fits of the peaks computed from Eq. (fTU]) for 
different values of N are shown in Fig. [31 The overall 
linewidth Aw, in general, decreases with the increase of 
current (Fig. [3^) , except in the region near the genera- 
tion threshold at which a local maximum of the linewidth 
can appear due to the lineshape asymmetry. Depend- 
ing on the magnitude of the nonlinear frequency shift 
N, this maximum can be several times larger than the 
generation linewidth in both subcritical and supercriti- 
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